Topological Change of the Fermi Surface in Low Density Rashba Gases: Application 

to Superconductivity 
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Strong spin-orbit coupling can have a profound effect on the electronic structure in a metal or 
semiconductor, particularly for low electron concentrations. We show how, for small values of the 
Fermi energy compared to the spin-orbit splitting of Rashba type, a topological change of the Fermi 
surface leads to an effective reduction of the dimensionality in the electronic density of states. We 
investigate its consequences on the onset of the superconducting instability. We show, by solving the 
Eliashberg equations for the critical temperature as a function of spin-orbit coupling and electron 
density, that the superconducting critical temperature is significantly tuned in this regime by the 
spin-orbit coupling. We suggest that materials with strong spin-orbit coupling are good candidates 
for enhanced superconductivity. 

PACS numbers: 72.25.-b, 72.10.-d, 72.20.Dp 
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Spin-orbit (SO) coupling arising from the lack of inver- 
sion symmetry plays a leading role in the field of spintron- 
ics One of the main goals in this field of research is the 
possibility of tuning the electron spin properties (trans- 
port, coherence, relaxation, etc.) by means of electrical 
fields [2J. With this aim in mind, different properties have 
been investigated, such as spin relaxation [3J], magneto- 
conductance Q and spin- Hall currents [f|. As a general 
rule, the spin-orbit coupling is assumed to be quite small 
with respect to the other relevant energy scales, in par- 
ticular with respect to the electronic dispersion, so that 
the infinite bandwidth limit is often employed. While 
this assumption is indeed rather reasonable in most of 
the cases, the natural aim of the current investigations is 
to search for new materials with stronger spin-orbit cou- 
plings, as for instance in HgTe quantum wells jf|, or the 
surface states of metals and semimetals 0, H|. For this 
reason, experimental evidence of a Rashba SO coupling 
with energy Eq (to be defined below) as large as ~ 220 
meV in bismuth/silver alloys @, or with E ~ 30 — 200 
meV in non-centrosymmetric superconductors CePtsSi 
HI [IH, Li 2 Pd 3 B, and Li 2 Pt 3 B [H [l3j], is certainly an 
important step towards the investigation of new materi- 
als with large SO coupling. Needless to say, the existence 
of such materials compels us to carry out a more thor- 
ough investigation of the properties of SO systems when 
Eq is no longer necessarily the smallest energy scale in 
the problem. 

The possibility of having novel interesting features in 
low density systems with Fermi energy Ep of the same 
order or lower than the SO energy Eq, in particular, has 
not been sufficiently investigated to date, in our opinion, 
and only few studies have been devoted to this problem. 
In Ref. [3] for instance, it was shown that the vanishing 



of the spin-Hall current in the low density limit n — > 
of Rashba disordered systems is not related to the van- 
ishing of the vertex function (which applies only in the 
strict E-p I Eq — * oo limit) but rather to the cancellation 
between on- Fermi surface and off-Fermi surface contribu- 
tions. Another interesting effect was also pointed out in 
Ref. 15j: there the spin relaxation time t s for E-p <C Eq 



was shown to be proportional to the electron scattering 
time r, in contrast with the standard Dyakonov-Perel 
behavior, where r s scales as 1/r [3]. On the other 
hand, in spite of the growing interest in the properties 
of non-centrosymmetric superconductors with strong SO 
coupling 13, [l8| , no specific investigation to explore the 
regime where Ep/Eo < 1 has been pursued, to our knowl- 
edge. 

The aim of this Letter is to explore in detail a fun- 
damental feature arising in SO Rashba systems, namely 
the topological change of the Fermi surface induced by 
the strong SO interaction in the low density regime 
Ep/Eo < 1- We show that in this situation the en- 
hanced phase space available for the electronic excita- 
tions gives rise to a SO-induced change of the electronic 
density of states which can be described in terms of an 
effective reduced dimensionality. We discuss the conse- 
quences of this scenario on the superconducting instabil- 
ity criterion for both two- and three-dimensional Rashba 
systems. We show that, in contrast with the high den- 
sity case Ep/E > 1, the SO coupling in the Ep/E < 1 
regime systematically enhances the superconducting crit- 
ical temperature T c , providing evidence that the lack of 
inversion symmetry can be remarkably beneficial for su- 
perconducting pairing. 

We begin our analysis by considering the Rashba model 
which describes the linear coupling of conduction 
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electrons with a spin-orbit potential of the form 

), (i) 

where <r x , <r y are Pauli matrices and 7 is the Rashba 
coupling constant. For two-dimensional systems such as 
asymmetric quantum wells and surface states of metals 
and semimetals the SO coupling in Eq. {T]) arises from 
the asymmetric confining potential, while in bulk three- 
dimensional compounds Eq. {T]) originates from the lack 
of reflection symmetry with respect to the z-direction, 
as in CePtaSi. The SO coupling is reflected in an en- 
ergy splitting of the two helicity bands. Assuming a 
parabolic band for 7 = 0, the resulting dispersion of 
the electronic excitations, for two-dimensional (2D) and 
three-dimensional (3D) cases, reduces to 

Ef(k) = ^(k±k„f, (2) 
E |B W . £ ( » ±W . + 0, (3) 

where k = |k| is the modulus of the xy in-plane mo- 
mentum k = (k x ,k y ), k z is the wavevector along the 
^-direction, and m* is the effective electron mass. In 
addition fco = m* r )/T% 2 represents here the character- 
istic Rashba momentum. The dispersion for the 2D 
case is shown in Fig. []Ji, which can easily be general- 
ized in the 3D case by taking into account the k z dis- 
persion. The two horizontal dashed lines correspond to 
the Fermi level for high density and low density regimes 
defined as E F > E and E F < Eq respectively, where 
Eq = h 2 kQ/2m* is the energy of the k = point with 
respect to the bottom band edge at k = k (see Fig. [1^). 

Density of states - Several studies in the literature have 
focused on the high density regime, E F Eq, where 
the two Fermi surfaces belong to different helicity bands, 
and where the Fermi volume V F is given by the area of 
two concentric Fermi circles Vf = 7rfcp . + irkp _ with 

/jf,± = \J 2m* I Ti 2 (V-E-f T VEq)- In this case the elec- 
tronic density of states (DOS) at the Fermi level is given 
by 

where the Fermi velocity |^f.±| = \J2E F jm* is indepen- 
dent of the helicity number s = ± and the Fermi sur- 
faces are <Sf.± — 2nkp,±. Hence the total DOS in the 
E F > £0 regime N 2D (E F ) = m*/(nh 2 ) is identical to 
the one in the absence of spin-orbit coupling. A similar 
result applies for the 3D case where, from Eq. ©, the 
corresponding DOS can be obtained as 

N 3B (E F ) = J ^N 2D (E F - h 2 k 2 z /2m*). (5) 




FIG. 1: Panel (a): electronic dispersion in the presence of 
SO coupling: for Ep > Eo (high density regime) the two 
Fermi surfaces belong to different helicity bands, while for 
-Bf < -Eo (low density regime) the Fermi surface exists only 
on the E^ band. Also shown here is the Rashba energy Eo 
corresponding to the lowest interband excitation energy for 
an electron lying at the bottom of the lower band. Panels (b) 
and (c): Fermi surface and density of states, respectively, in 
the low density regime for the two-dimensional case. Panels 
(d) and (e): Fermi surface and density of states in the low 
density regime for the three-dimensional case. 

In the high density regime E F > Eq we get N 3I) (E F ) = 
a{\/E F - Eq + \/E~q~ arctan[- v /i?o/(£ , F - ^o)]}, where a = 
y/2m* 3/2 /(ir 2 h 3 ), which reduces to N 3B (E F ) ~ a^/Ep' in 
the E f /Eq ^S> 1 limit. This again is the result one obtains 
in the absence of SO coupling. 

Let us now consider the E F < Eq regime. In this 
case the Fermi level intersects only the lower E_(k) 
band and the topology of the Fermi surfaces dras- 
tically changes. In the 2D case for instance only 
the annulus that lies be tween two Fermi circles of 

radii fop, 2 = "J 2m* I a 2 {\[Eq~ + y/E-p) and fcF.i = 

\j2m*/h 2 {\/Eq — t/E f ) , belonging to the same helicity 
band, is filled (Fig. QJ>), and the inner Fermi surface is 
inwards oriented. We can still employ Eq. ([4]) by sum- 
ming over the two Fermi surface indexes s = 1,2. Using 
|t>F, s | = \/2E F Jm* and S FyS = 2irk F . s we get 

^ 2D (£f) = ^#, (6) 
nn V -&F 

which is valid as long as E F < Eq (Fig. [TJ;). Most pe- 
culiar is the square-root divergence for E F — > that is 
reminiscent of one-dimensional behavior. We relate such 
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a feature to the non-vanishing in the low density limit 
of the Fermi surface which remains finite, Sp <a oc \/~Eq, 
while at the same time, the Fermi velocity vanishes as 
\fE-p. The behavior of Eq. © has to be compared with 
the 7 = case where also the Fermi surface shrinks as 
\fEp and the electron DOS has a non divergent step-like 
behavior in the Ep — > limit [20I ]. 

A similar reduction of effective dimensionality in the 
electron DOS appears also for the 3D systems in the 
Ep < Eq regime. In this case the Fermi surface has a 
torus-like topology as shown in Fig. [TJi, with major ra- 
dius fco = *J 2m*E /h 2 and minor radius J 2m*Ep/h 2 . 
Applying once more Eq. (J3J) we get 



N 3D (E F ) = ^V^o 



(7) 



for Ef < Eq. We see therefore that the SO coupling- 
changes qualitatively the low density behavior of the 
DOS providing a finite step-like behavior (Fig. [T^) in 
contrast with a standard 3D electron gas whose DOS 
vanishes as y/Ep. 

Cooper instability - The above described reduction of 
effective dimensionality sheds a new light on the possible 
existence of a superconducting phase in the low density 
regime of systems with no inversion symmetry. To illus- 
trate this point let us consider the classical problem [2l[ 
of a Fermi surface instability towards the formation of a 
Cooper pair: 



1 = V 



1 



2£ + A' 



(8) 



where A > is the binding energy of the bound pair 
state, V is the strength of the interaction which we con- 
sider here for simplicity in the s-wave channel, and where 
we have introduced a standard BCS cut-off ujq related 
to the energy of the underlying bosonic mediator. Since 
the superconducting Cooper pairing is essentially a Fermi 
surface instability, the strength of the bound state and its 
existence itself is intimately related to the phase space of 
the available electronic excitations. For instance, as well 
known, in the low density limit of 3D systems, where 
N (£) = ay/%, Eq. ijSJ) predicts a finite critical coupling 
V c = l/(a^/UJo) below which no bound state exists. 

This result changes drastically for finite SO couplings 
where, as seen above, the electron DOS behaves now as 
an effective 2D system with a constant value in the low 
energy regime. Using Eq. in Eq. © we get: 



A 3D ~ 2u> exp 



naVJE, 



(9) 



which explicitly shows that, contrary to the usual 3D 
case, the Cooper pair instability exists no matter how 
weak V is. Furthermore Eq. |9j) predicts that the pair 
energy A has an exponential dependence on the SO cou- 
pling, as is normally the case (here as well) for the at- 
tractive potential. 



A similar change of the character of the Cooper insta- 
bility occurs also in the 2D case. Indeed, in the absence 
of SO coupling, one would get the standard BCS result 
A = 2ujq exp(— 2nh 2 /m*V). On the other hand, due to 
the strong one-dimensional- like divergence of the electron 
DOS, Eq. ©, the binding energy for finite SO coupling 
reads now 



1 

A 2D = - 



i*V 



(10) 



where the bosonic energy ujq is no longer present and 
the relevant energy scale is provided by Eq. Note also 
the quadratic dependence of the binding energy A with 
respect to the coupling strength V, and the complete 
absence of an isotope effect for the phonon- mediator case. 

Superconducting critical temperature - The above dis- 
cussion of the single Cooper pair problem will be now 
a guide to the following investigation of the supercon- 
ducting transition for finite (low) densities in fully in- 
teracting systems. To this end we consider a Rashba- 
Holstein model where the SO coupled electrons interact 
with dispersionless bosons with ener gy u jq through an s- 
wave coupling with matrix element g [22| . The supercon- 
ducting properties, and in particular the critical temper- 
ature T c , are evaluated within the Eliashberg framework 
properly generalized in the presence of SO coupling, with 
general values of Ef/Eq, including the low density case 
just discussed. Note that due to the lack of inversion 
symmetry, a different superconducting phase with mixed 
even/odd order parameter and singlet /triplet symmetry 
can in principle be included [l8j]. We neglect here for 
simplicity this issue and focus only on the s-wave singlet 
channel. 

In Fig. [2^,c we show the superconducting critical tem- 
perature T c as a function of the SO Rashba energy Eq 
for different electron densities n. In the figure the lower 
density values, ) 



10 13 cm" 2 and n = 10 2 



respond to Fermi energies Ep ~ 24 meV and Ep ~ 46 
meV for the free electron gas in 2D and 3D respectively. 
In addition, for a practical purpose one needs to intro- 
duce a finite bandwidth cut-off E c , which is physically 
provided by the size of the Brillouin zone k c . We set 
E c = 2000(430) meV which give k c ~ 0.72(0.33) A" 1 
respectively for the 2D and 3D case, We get thus an 
energy-dimensional electronic DOS per unit cell which, 
for the density values reported above in the absence of 
SO coupling , is N 2 d(E f = 24meV) ~ 5 • 10" 4 meV" 1 
and N 3D (E F = 46meV) ~ 12 • 10~ 4 meV -1 . For all cases 
ujq has been fixed at ujq = 20 meV and g = 5wo- With 
these values we obtain dimensionless coupling constants 
A = 2g 2 N(Ep)/ujQ respectively A 2D (£ F = 24meV) ~ 0.5 
and A3d(^f = 46meV) ~ 0.6. Fig. [2ji,c shows a signifi- 
cant increase of T c as a function of Eq, in particular for 
low densities where a small Eq is sufficient to enter into 
the Ep < Eq regime. This holds true for both 2D and 
3D systems, and the enhancement of T c can be as high as 
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FIG. 2: Panels (a) and (b): Superconducting critical tem- 
perature T c as function of the Rashba spin-orbit energy Eq 
and of the electron density n, respectively, for the 2D case. 
Panels (c) and (d): same quantities for the 3D system. 



300 % with respect to the Eq — > limit. Also interesting 
is the study of the superconducting critical temperature 
as a function of the electron density n, as reported in 
Fig. [2b, d. In these panels it is evident how the T c vs. 
n behavior reflects the effective dimensionality reduction 
of the underlying non interacting DOS. For example for 
the 2D case T c is nearly constant as a function of electron 
density for Eq = 0; this is consistent with the constant 
electron DOS. In contrast, for finite Eq the presence of 
strong peaks in T c vs. n reflects the ID-like singularity of 
the DOS. Note however that the retarded electron-boson 
interaction gives rise to dynamical one-particle renormal- 
ization effects, automatically taken into account in the 
Eliashberg equations, which smear the singularity of the 
bare DOS. 

Similarly, in the Eq — > limit of the 3D case T c 
drops as the density n is reduced due to vanishing of 
Nsd(E-p) oc \ZE-p. On the other hand the 2D charac- 
ter of the DOS with Eq ^ gives rise to an almost flat 
dependence of T c for sufficiently low n, with a critical 
temperature tuned by the Rashba energy Eq. Both the 
2D and 3D case thus show that the lack of inversion sym- 
metry not only affects the character of the order param- 
eter, as discussed in several works 
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SO interaction. 

Let us discuss now the relevance of our results in the 
context of real materials. Concerning the 2D case, for 
instance, surface states and low dimensional heterostruc- 
tures could be natural candidates for the search of en- 
hanced superconductivity. In particular the issue of sur- 
face superconductivity [23j has recently been brought to 
attention due to its relevance for systems like alkali-doped 
WO 3 [3] where, for sufficiently low concentrations of 
the alkali atoms, evidence of superconductivity confined 
to the surface has been provided; this is precisely where 
strong SO coupling arises from the confinement of the 
surface potential. Interesting perspectives are also given 
by the non-centrosymmetric superconductors, such as 
CePtsSi, Li2Pd3B, and I^PtaB, where a strong Rashba 
energy arises from the lack of inversion symmetry in the 
bulk crystal. In the I^PdsB and L^PtaB compounds, 
in particular, a large SO coupling is accompanied by a 
strong electron-phonon interaction (25|. In this case, of 
course, as in the heavy fermion CePtaSi case, a proper 
generalization of the present results to the case of non 
parabolic bands is needed. 

In summary, we have examined the impact of a strong 
Rashba spin-orbit interaction on superconductivity in a 
weakly coupled electron-boson system. The primary re- 
sult is an effective reduction in dimensionality of the 
electronic density of states. This accomplishes what the 
Fermi sea did for the Cooper pair problem in 3 dimen- 
sions — it increased phase space at the Fermi level sig- 
nificantly to allow binding for arbitrarily weak interac- 
tions. We then performed full Eliashberg calculations 
of the critical temperature for a wide range of parame- 
ters; these illustrate that the enhanced density of states 
has a significant impact on the superconducting critical 
temperature. We suggest a search for higher critical tem- 
peratures in materials with large spin-orbit coupling. In 
systems where electron density can be varied one should 
be able to test some of the trends reported here. 

This work was supported in part by the Natural 
Sciences and Engineering Research Council of Canada 
(NSERC), by ICORE (Alberta), by the Canadian Insti- 
tute for Advanced Research (CIAR), and by the Univer- 
sity of Geneva. 



principle can also lead to a substantial enhancement of 
the superconducting pairing in the low density regime. 
As explained above, this phenomenon is triggered by the 
topological change of the Fermi surface due to the strong 
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